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I  100)  at  a  period  of  100  s  than  it  is  a  1  s  1000). 
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explained  H^y  fluid  flow  in  a  network  of  cracks  in  brittle  rock.  Increasing  Q  with  depth 
to  10-15  km  is  caused  by  the  increased  closing  of  those  cracks  due  to  pressure  and  their 
enhanced  healing  with  increasing  temperature.  Viscous  flow  at  greater  depths  contributes 
further  to  the  dissipation  of  cracks  and  to  further  increases  in  Q  .  Decreasing 
values  at  greater  depths  can  be  explained  as  the  result  of  increasing  temperature, 
increasing  content  of  partial  melt,  enhanced  dislocation  motion,  or  some  combination  of 
these  effects. 

An  Lg  coda  Q  map  has  been  derived  for  nearly  all  of  Eurasia.  Q  values  range  from 
1000  or  more  in  stable  region  to  200  or  less  in  regions  having  undergone  recent 
tectonic  activity.  Q  values  in  different  shield  regions  differ,  the  maximum  values 
being  lower  in  the  Siberian  platform  than  those  in  the  East  European  platform. 
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TECHNICAL  SUMMARY 

A  new  inversion  method,  in  which  the  frequency  dependence  of  shear  wave 
internal  friction  is  allowed  to  vary  with  depth,  was  developed  and 

applied  to  selected  Rayleigh  wave  attenuation  data  in  the  Basin  and  Range 
Province  of  the  southwestern  United  States,  Models  were  further  constrained 
to  satisfy  observed  values  of  Q  and  the  frequency  dependence  of  Lg  waves  in 
that  region.  Many  models  can  explain  those  two  data  sets  within  their  uncer¬ 
tainties,  but  at  1  Hz  most  have  low  values  of  (50-80)  in  the  upper  8  km  of 
the  crust,  rapidly  increasing  values  to  about  1000  at  mid-crustal  depths,  and 
decreasing  values  at  greater  depths.  Models  which  include  a  layer  of  higher 
Q  values  (80-150)  in  the  upper  few  km  of  the  crust,  overlying  a  region  of  low 
Q  values,  cannot  be  precluded  by  the  data  of  this  study.  Models  for  which 
the  frequency  dependence  is  low  (0.0-0. 1)  in  the  upper  crust  best  explain 
the  available  data.  In  the  lower  crust  that  frequency  dependence  may  be 
substantial  and  leads  to  models  in  which  may  be  an  order  of  magnitude 
lower  (^  100)  at  a  period  of  100  s  than  it  is  a  1  s  (^  1000). 

The  low  values,  and  their  constancy  with  frequency,  in  the  upper¬ 
most  crust  can  be  explained  by  fluid  flow  in  a  network  of  cracks  in  brittle 
rock.  Increasing  with  depth  to  10-15  km  is  caused  by  the  increased 
closing  of  those  cracks  due  to  pressure  and  their  enhanced  annealing  with 
increasing  temperature.  Viscous  flow  at  greater  depths  contributes  further 
to  the  dissipation  of  cracks  and  to  further  increases  in  .  Decreasing 
values  at  greater  depths  can  be  explained  as  the  result  of  increasing  tempera¬ 
ture,  increasing  content  of  partial  melt,  enhanced  dislocation  motion,  or  some 
combination  of  these  effects. 
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An  Lg  coda  Q  map  has  been  derived  for  nearly  all  of  Eurasia.  Q 
values  range  from  1000  or  more  in  stable  regions  to  200  or  less  in  regions 
having  undergone  recent  tectonic  activity.  Q  values  in  different  shield 
regions  differ,  the  maximum  values  being  lower  in  the  Siberian  platform 
than  those  in  the  East  European  platform. 
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1  INTRODUCTION 

Records  of  surface  waves  at  regional  distances  in  continental  regions  typically 
consist  of  the  short-period  Lg  phase,  along  with  its  coda,  and  longer  period 
fundamental-mode  waves.  Modern  broad-band  digital  seismic  stations  record 
these  multiphase  surface  waves  on  the  same  data  channel,  even  though  their 
frequency  contents  are  quite  different.  We  presented,  in  earlier  papers, 
results  of  our  investigations  of  the  dynamic  properties  of  Lg  and  Lg  coda  (Xie 
and  Mitchell,  1990)  and  of  fundamental-mode  surface  waves  (Mitchell  et  al., 
1992)  and  showed  how  they  are  affected  by  laterally  complex  structure  in  a 
region  of  the  southwestern  United  States  which  includes  the  Basin  and  Range 
and  adjacent  tectonic  provinces.  These  papers  will  hereafter  be  referred  to  as 
Part  I  and  Part  II,  respectively. 

When  propagating  through  a  laterally  complex  structure  the  dynamic 
properties  of  Lg  and  fundamental-mode  surface  waves  are  affected  by  various 
types  of  scattering  produced  by  that  complexity,  as  weU  as  by  intrinsic 
absorption  of  the  medium.  Ideally,  we  would  like  to  isolate  data  which  are 
affected  solely,  or  at  least  predominantly,  by  intrinsic  absorption,  or  anelasti- 
city.  The  Q  values  which  quantify  that  absorption  are  determined  by  such 
factors  as  temperature,  and  the  amounts  of  fluid  content,  partial  melt,  or 
solid  state  dislocation  motion  in  the  crust  or  upper  mantle.  Knowledge  of 
intrinsic  Q  values  within  the  crust  and  upper  mantle  may  therefore  contribute 
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to  our  understanding  of  the  tectonic  evolution  of  regions  where  those  values 
have  been  determined. 

The  acquisition  of  data  sets  which  are  suitable  to  invert  for  intrinsic 
anelasticity  or  Q  requires  great  care.  The  sensitivity  of  Q  measurements  to 
particular  kinds  of  lateral  complexity  can  be  quite  different  for  different 
phases.  It  is  important,  therefore,  when  using  these  phases  to  study  crustal 
Q  structure,  to  understand  possible  biases  which  can  affect  them,  and  to 
either  discard,  or  apply  corrections  to,  data  which  are  adversely  affected  by 
those  biases. 

Part  I  analyzed  short-period  Lg  and  Lg  coda  recorded  by  stations  of  the 
Lawrence  Livermore  National  Laboratory  (LLNL)  network  and  resulted  in 
highly  consistent  Lg  Q  (QLg)  Lg  coda  Q  (Q^)  values  for  a  propagation 
path  across  the  Basin  and  Range  province.  This  was  achieved  using  a  care¬ 
fully  designed  inverse  scheme  and  a  data  selection  procedure  in  which  we 
reduced  the  effects  of  focusing/defocusing  and  non-flat  site  response  on  the 
Lg  phase,  and  excluded  those  coda  data  which  were  significantly  contam¬ 
inated  by  scattering  or  by  reflections  from  the  coast  line  outside  the  LLNL 
network  (Kennett  et  al,  1990).  Part  n  studied  intermediate-period 
fundamental-mode  surface  wave  propagation  and  showed  that  particle 
motion,  and  hence  surface-wave  amplitudes,  were  often  adversely  affected  by 
major  tectonic  features  such  as  coast  lines  or  the  Sierra  Nevada,  especially  if 
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the  stirface  waves  impinged  obliquely  on  those  features.  One  path  (ELK- 
MNV),  however,  which  is  sub-perpendicular  to  the  Sierra  Nevada,  was  found 
to  provide  unbiased  data. 

The  present  paper  utilizes  data  from  Part  I  and  Part  11,  which  were 
found  to  be  unbiased,  to  invert  for  crustal  anelasticity  structure  for  the  Basin 
and  Range.  Implicit  in  our  inversion  procedure  is  the  assumption  that  Q 
measurements  for  both  types  of  data  are  due  to  intrinsic  shear-wave  internal 
friction  (  Q~^  )  rather  than  to  scattering.  That  this  assumption  is  satisfied  for 
long-period  Rayleigh  waves  is  apparent  from  the  observation  in  Part  11  that 
retrograde-elliptical  particle  motion  characterizes  those  waves  along  the 
selected  path  of  study.  QLg  and  its  frequency  dependence  along  the  selected 
path  are  also  assumed  to  have  been  produced  by  intrinsic  anelasticity  effects. 
This  conclusion  is  based  upon  the  equivalence  of  QLg  and  Qj.  in  Part  I  and  on 
recent  studies  (e.g.  Frankel  and  Wennerberg,  1987)  which  have  found  that 
coda  Q  is  representative  of  intrinsic  Q.  Further  support  of  that  assumption  is 
provided  by  the  work  of  Jin  and  Aki  (1988)  who  found  that  coda  Q  in  a  low- 
Q  region  of  China  is  likely  to  be  measuring  intrinsic  Q  and  by  Kennett  (1990) 
who  demonstrated  theoretically  that  scattering  contributes  considerably  less 
to  apparent  QLg  for  a  low  Q  environment  than  it  does  for  a  high-Q  environ¬ 
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This  paper  attempts  to  jointly  utilize  intermediate-period  (6-33  s) 
fundamental-mode  Rayleigh  waves  and  QLg  at  1  Hz  to  obtain  a  model  of 
shear  wave  internal  friction  for  the  crust  and  uppermost  mantle  of  the  Basin 
and  Range  province.  We  will  utilize  a  new  inversion  method  in  which  the 
frequency  dependence  of  Q^,  is  allowed  to  vary  with  depth  and  we  will 
assume  that  both  Rayleigh  wave  attenuation  and  QLg  are  dominated  by 
intrinsic  anelasticity  rather  than  by  scattering  effects. 

2  PREVIOUS  CRUSTAL  Q  STUDIES  IN  THE  BASIN  AND  RANGE 

Previous  studies  of  Q  in  the  Basin  and  Range  province  using  various  wave 
types  have  aU  indicated  that  crustal  Q  values  there  are  lower  than  in  stable 
regions  of  North  America.  Sutton  et  al.  (1967)  found  the  first  indication  of 
low  Q  from  the  relative  fall-offs  of  Pg  and  Lg  energies  for  events  in  the  Basin 
and  Range  and  in  the  central  United  States.  Mitchell  (1975)  obtained  a  model 
of  crustal  shear  wave  Q  (Q^,)  for  the  western  United  States  and  found  that 
differences  in  Rayleigh  wave  attenuation  between  the  eastern  and  western 
United  States  could  be  explained  by  lower  values  of  Q^.  in  the  upper  crust  in 
the  western  United  States.  Cheng  and  Mitchell  (1981),  using  a  single-station 
multi-mode  approach,  found  values  in  the  upper  crust  of  the  Basin  and 
Range  to  be  about  85,  as  compared  to  about  275  in  the  central  United  States. 
Patton  and  Taylor  (1984)  inverted  Rayleigh  wave  attenuation  data  obtained 
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over  a  broad  portion  of  the  Basin  and  Range  and  presented  three  models 
in  which  values  in  the  upper  crust  ranged  between  about  110  and  170  and 
values  in  the  lower  crust  which  were  substantially  lower  than  those  obtained 
by  Mitchell  (1975)  and  Cheng  and  Mitchell  (1981).  Lin  (1989)  inverted 
attenuation  data  obtained  for  a  two-station  path  in  the  western  Basin  and 
Range  and  obtained  Q^.  values  between  about  40  and  60  in  the  upper  crust 
and  high  values  (1000  or  more)  at  mid-crustal  depths.  A  similar  model 
obtained  by  Mitchell  (1991)  explains  both  fundamental  Rayleigh  mode 
attenuation  and  observed  Lg  Q  values. 

A  finding  common  to  aU  of  the  crustal  Q^,  determinations  in  the  Basin 
and  Range  province  is  that  the  upper  crust  has  substantially  lower  Q^,  values 
than  those  found  in  more  stable  regions.  This  low-Q  layer  provides  an  expla¬ 
nation  for  the  high  relative  attenuation  of  Lg  reported  for  that  region  by  Sut¬ 
ton  et  al.  (1967)  and  is  supported  by  subsequent  determinations  of  Q  for  the 
Lg  phase  (QLg)  arid  Lg  coda  (Q^)  which  are  lower  than  for  stable  regions  (e.g. 
Singh  and  Herrmann,  1983;  Peseckis  and  Pomeroy,  1984;  Chavez  and  Pri¬ 
estly,  1986;  Nuttii,  1986). 

3  FREQUENCY  DEPENDENCES  OF  QLg  AND 

An  expression  for  the  internal  friction  of  Rayleigh  waves  (Qr^)  was  given  by 
Anderson  et  al.  (1965).  Using  that  expression,  Qr^  can  be  obtained  for  any 
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Rayleigh  mode  if  the  appropriate  velodties,  partial  derivatives  of  Rayleigh 
wave  phase  velocities  with  respect  to  shear-  and  compressional-wave  veloci¬ 
ties,  and  intrinsic  values  are  used.  Mitchell  (1980)  extended  that  expres¬ 
sion  to  incorporate  Q^,  values  which  vary  with  frequency.  It  is  sometimes 
assumed  that  this  frequency  dependence  of  also  pertains  to  QLg  (e.g. 
Campillo  ei  al.,  1985;  Herraiz  and  Espinosa  1987;  Nuttli,  1986).  Mitchell  (1991) 
pointed  out,  however,  that  the  frequency  dependence  of  Qr~^  can  be  caused 
by  layering  of  the  Q^,  structure  in  the  crust  as  weU  as  by  intrinsic  frequency 
dependence  of  Q^.  He  further  showed  that  the  frequency  dependence  of  Lg 
observed  in  the  Basin  and  Range  province  can  be  entirely  explained  by  a  lay¬ 
ered  model  in  which  is  independent  of  frequency  and  increases  from  low 
to  high  values  at  mid-crustal  depths. 

Frequency  dependences  of  Q^.  and  Qi^g  are  therefore  fundamentally  dif¬ 
ferent  parameters.  In  the  sections  that  follow  ^  will  refer  to  intrinsic  fre¬ 
quency  dependence  of  whereas  ti  refers  to  frequency  dependence  of  QLg 
which  is  produced  by  the  combined  effects  of  intrinsic  i  and  changing  values 
of  with  increasing  depth  in  the  crust. 

4  DATA 

Data  used  to  invert  for  crustal  anelasticity  in  this  study  consist  of 
fundamental-mode  Rayleigh  wave  attenuation  coefficient  values  from  Part  n 


-8- 


for  the  path  between  stations  ELK  and  MNV  (Figure  1).  That  path  is  sub¬ 
perpendicular  to  the  major  structural  trends  in  the  region  and  particle  motion 
plots  show  that  those  values  are  not  seriously  affected  by  lateral  refraction  or 
multipathing  of  the  Rayleigh  waves.  Those  data  were  discussed  in  detail  in 
Partn. 

Further  constraints  on  models  obtained  from  the  inversion  of 
fundamental-mode  data  are  provided  by  the  Lg  Q  values  and  their  frequency 
dependence  obtained  in  Part  I.  After  removing  data  for  which  Lg  and  its 
coda  were  adversely  affected  by  portions  of  oceanic  path,  a  very  stable  esti¬ 
mate  was  found  for  QLg  and  its  frequency  dependence  in  the  Basin  and 
Range  province.  Using  the  relation  QLg(f)  ®  Qo^/  where  Qq  is  Qi^g  at  1  hz 
and  T|  is  a  frequency  dependence  exponent,  it  was  found  in  that  QLg(0  “  (267 
±  Since  these  values  for  and  iq  are  very  similar  to  those 

found  for  coda  Q  along  the  same  path,  we  inferred  in  Part  I  that  they  were 
produced  by  intrinsic  anelasticity,  rather  than  by  scattering.  Any  model 
obtained  from  the  inversion  of  fundamental-mode  Rayleigh  wave  data  should 
therefore  predict  a  Qq  value  of  about  270  and  a  Qi^g  frequency  dependence  of 
about  0.4. 

5  INVERSION  PROCEDURE  FOR 

The  discussions  in  the  previous  section  and  in  Parts  I  and  n  suggest  that  the 
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attenuation  of  both  fundamental-mode  Rayleigh  waves  and  1-Hz  Lg  waves  in 
the  crust  of  the  Basin  and  Range  province  are  dominated  by  intrinsic  absorp¬ 
tion.  The  same  model  of  shear  wave  internal  friction  for  that  region 

should  therefore  explain  both  types  of  data. 

Mitchell  (1980)  modified  the  equations  of  Anderson  et  al  (1965)  to  allow 
to  vary  with  frequency,  f,  but  restricted  that  frequency  variation  to  be 
constant  with  depth.  In  the  present  study  we  allow  that  frequency  depen¬ 
dence  to  vary  with  depth,  so  it  can  be  written 


Q,i  =  Qo/'  (1) 

where  the  subscript  1  is  a  layer  index  and  Qqi  is  the  Q^j,  value  at  1  Hz  for  layer 

1.  Toksoz  et  al.  (1990)  have  previously  allowed  the  frequency  dependence  of 
shear  wave  Q  to  vary  with  depth  in  a  study  in  which  they  compared  theoreti¬ 
cal  and  observed  seismograms  at  higher  frequencies  to  infer  an  atttenuation 
model  for  the  crust  in  a  high-Q  region. 

We  invert  attenuation  coefficient  data  (7r)  for  using 
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where  the  subscripts  R,  ot  and  3  identify  Rayleigh,  compressional,  and  shear 
waves,  respectively,  o)  is  angular  frequency,  and  Cr  is  Rayleigh  wave  phase 
velocity.  When  co,  a,  and  3  appear  as  subscripts,  they  indicate  that  those 
quantities  are  held  constant  when  the  partial  derivatives  are  computed.  The 
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factor  —  occurs  in  the  right-hand  term  because  we  assume  that  compressional 


wave  Q  is  twice  as  large  as  for  all  inversions  of  this  study.  Previous  stu¬ 
dies,  such  as  Mitchell  (1980),  have  indicated  that  inversion  results  are  not 
strongly  affected  by  the  assumed  value  of  compressional  wave  Q.  Values  for 
tti  and  Pi  are  taken  from  the  Basin  and  Range  model  of  Lin  (1989)  and  values 
3C  3C 

of - , - ,  and  Cr  are  computed  for  that  model. 

3ai  3  Pi 


The  inversion  process  proceeds  as  follows: 

1.  Assume  a  depth  distribution  of  ^i  and  invert  fundamental-mode 
Rayleigh-wave  attenuation  coefficients  to  obtain  a  model  for  the 
Basin  and  Range. 

2.  Calculate  the  fundamental-mode  Rayleigh  wave  attenuation  coefficient 
values  which  are  predicted  by  that  model  and  compare  with  observed 
values.  If  they  agree  within  the  uncertainties  of  the  data,  go  on  to  3;  if 
not,  try  a  new  distribution  of  Ci- 

3.  Compute  synthetic  Lg  seismograms  at  several  distances  from  a  seismic 
source  using  the  velocity  model  of  Lin  (1989)  and  the  model  result¬ 
ing  from  the  inversion  of  the  fundamental-mode  attenuation  coefficient 
data.  An  example  of  several  such  S5mthetics  at  distances  between  400 
and  1200  km  appears  in  Figure  2.  In  the  inversions  of  the  present  study, 
we  computed  20  seismograms  evenly  spaced  at  distances  between  300 
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and  1250  km. 

4.  Apply  the  stacked  spectral  ratio  (SSR)  method  of  Xie  and  Nuttli  (1988)  to 
the  set  of  S5mthetic  Lg  seismograms  to  obtain  values  of  Qq  and  tj 
predicted  by  the  derived  Qj];;'  model,  and  compare  those  values  wih 
observed  values  of  Qq  and  t|  in  the  Basin  and  Range. 

This  process  can  be  repeated  from  several  different  distributions  of  Cj  to 
obtain  those  models  which  best  explain  both  the  fundamental-mode  Rayleigh 
wave  attenuation  data  and  the  values  of  Qq  and  t|  for  Lg  waves. 

Our  inversion  procedure  assumes  a  layered  structure  in  which  the 
number  of  layers  exceeds  the  number  of  observed  attenuation  values.  This 
makes  the  inversion  of  Rayleigh  wave  attenuation  underdetermined  and 
therefore  non-unique.  Non-uniqueness  may  be  partly  responsible  for  the 
differences  between  previously  obtained  models  of  Q^,  (Patton  and  Taylor, 
1984;  Lin,  1989).  The  degree  of  non-uniqueness  in  the  present  study  is 
higher  than  those  in  previous  studies  since,  by  allowing  the  frequency  depen¬ 
dence  of  to  vary  with  depth,  we  have  introduced  another  set  of  parame¬ 
ters  to  the  inversion  process.  Because  of  this  non-uniqueness  we  have  stu¬ 
died  three  classes  of  models  to  see  how  well  they  explain  fundamental-mode 
Rayleigh  wave  attenuation  and  the  Q  and  frequency  dependence  of  1-Hz  Lg 
waves.  The  three  classes  of  models  are  those  in  which  (1)  t,  is  uniform  with 
depth,  (2)  t,  increases  from  0.0  in  the  upper  crust  to  higher  values  at  depth. 
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and  (3)  £  assumes  positive  values  in  the  upper  crust  and  decreases  to  0.0  at 
depth. 

An  example  of  a  model  resulting  from  the  inversion  of  Rayleigh  wave 
attenuation  data  for  the  case  of  ^  =  0.0  at  all  depths  appears  in  Figure  3.  The 
effect  which  this  model  has  on  QLg  is  discussed  in  the  following  section.  The 
model  is  characterized  by  low  Q^,  values  (~  50)  in  the  uppermost  crust,  high 
values  (~  1000)  at  mid-crustal  depths,  and  decreasing  values  at  greater 
depths.  The  resolving  kernels  to  the  right  of  the  model  indicate  that  resolu¬ 
tion  is  best  in  the  upper  crust  and  degrades  with  increasing  depth.  At  any 
depth,  the  value  obtained  is  an  average  which  is  smeared  out  over  the  depth 
range  of  the  resolving  kernel.  Those  kernels  indicate  that  only  the  gross 
features  of  the  models  can  be  resolved  and  any  changes  which  occur  over 
small  depth  ranges  have  no  real  meaning.  The  standard  deviation  bars  asso¬ 
ciated  with  the  model  are  small  enough  so  that  we  can  attribute  significance 
to  the  major  features.  There  is,  of  course,  a  preset  parameter  which  allows  us 
to  trade  off  between  resolution  and  standard  deviation  in  the  inversion  pro¬ 
cess.  We  tried  several  settings  for  that  parameter  in  each  inversion  and 
selected  the  model  which  best  fit  the  Rayleigh  wave  attenuation  data  and,  at 
the  same  time,  did  not  lead  to  negative  values  or  wildly  fluctuating 
changes  between  layers.  Since  all  models  of  this  study  were  obtained  in  this 
way  the  widths  of  the  resolving  kernels  and  standard  deviations  for  all  of 
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them  are  similar  to  those  for  the  model  in  Figure  3. 

6  INVERSION  RESULTS 

Models  were  first  obtained  in  which  the  frequency  dependence,  5/  of  Qjt  was 
held  at  a  uniform  value  for  all  depths.  Three  of  those  models,  pertaining  to  a 
frequency  of  1  Hz,  appear  in  Figure  4  for  ^  =  0.0,  0.1,  and  0.2.  All  three 
models  adequately  explain  fundamental-mode  attenuation  data  (Figure  5), 
although  the  model  for  which  ^  =  0.2  would  predict  an  attenuation  coefficient 
value  which  is  slightly  outside  the  standard  deviation  bar  at  6  s.  Values  of 
Qo  and  t]  predicted  by  the  models  also  appear  in  the  figure.  The  model  in 
which  C  =  0.1  provides  the  best  agreement  with  the  observed  value  of  Qo 
(267)  but  predicts  t|  values  which  are  somewhat  higher  than  the  observed 
value  of  0.37.  The  model  in  which  i  =  0.0  provides  the  best  agreement  with 
the  observed  -ri  value  but  predicts  a  somewhat  low  Qo  value.  Figure  5  com¬ 
pares  values  of  Rayleigh  wave  attenuation  coefficients  predicted  by  models 
with  t,  values  of  0.0,  0.5,  and  1.0  with  observed  values.  While  all  models 
predict  attenuation  coefficient  values  which  lie  within  one  standard  deviation 
of  mean  observed  values  at  periods  greater  than  8  s,  only  the  curve  for  the 
model  in  which  ^  =  0.0  lies  within  the  standard  deviations  of  observed  values 
at  shorter  periods.  These  comparisons  suggest  that  ^  for  this  class  of  models 
must  be  about  0.2  or  less.  The  Qo  and  if]  values  exhibited  in  Figure  4  also 
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imply  that  for  models  with  a  constant  Q^,  frequency  dependence  with  depth, 
^  must  be  no  greater  than  about  0.2. 

The  next  group  of  models  which  we  studied  are  those  in  which  the  fre¬ 
quency  dependence  of  is  0.0  at  shallow  depths  and  assumes  a  higher 
value  deeper  in  the  crust.  It  was  possible  to  find  several  models  which 
predict  and  ir|  values  which  are  consistent  with  observed  values.  Figure  6 
shows  three  of  these  and  Figure  7  indicates  how  weU  they  agree  with 
observed  Rayleigh  wave  attenuation  coefficient  values.  A  model  in  which  ^  = 
0.0  at  depths  between  0  and  15  km  and  ^  =  0.5  at  greater  depths  produced 
and  Tj  values  which  are  in  excellent  agreement  with  observed  values.  When  I 
at  depth  is  increased  to  1.0  a  model  in  which  that  increase  occurs  at  a  depth 
of  18  km  also  provides  good  agreement  between  predicted  and  observed 
values  of  Qq  and  ti  for  Lg  waves.  When  the  ^  increase  is  only  to  0.2  the  best 
fitting  model  is  that  for  which  the  increase  occurs  at  a  depth  of  8  km.  That 
model  predicts  a  Qq  value  which  is  in  excellent  agreement  with  observed  Qq 
but  it  predicts  an  iq  value  which  is  somewhat  high.  The  models  in  which  i 
increases  to  0.2  and  to  0.5  both  predict  Rayleigh  wave  attenuation  coefficient 
values  which  fall  within,  or  are  very  close  to,  aU  observed  standard  devia¬ 
tions.  The  model  in  which  ^  increases  to  1.0  at  a  depth  of  18  km  predicts  a 
Rayleigh  wave  attenuation  coefficient  values  which  are  somewhat  low  at 
periods  of  6  and  8  s,  but  agree  with  observed  values  at  aU  longer  periods. 
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These  results  indicate  that  many  models  in  which  ^  is  small  near  the  surface 
and  increases  at  depth  in  the  crust  can  explain  both  fundamental-mode  Ray¬ 
leigh  wave  attenuation  coefficient  data  and  Qq  and  in  values  for  Lg  quite  well. 
The  degree  of  non-uniquess  in  these  inversions  is,  therefore,  fairly  large;  a 
range  of  ^  values  can  produce  similar  results.  In  addition  the  depth  at  which 
$  increases  from  zero  to  a  positive  value  is  not  well  resolved.  For  instance,  if 
^  increases  from  0.0  to  0.5,  the  depth  at  which  that  increase  occurs  can  vary 
between  about  10  and  about  20  km  and  still  satisfy  fundamental-mode 
attenuation  data,  as  well  as  Qq  and  in  values  within  their  uncertainties. 

Another  group  of  models  was  studied  in  which  £  assumed  a  positive 
value  at  shallow  depths  and  changed  to  zero  at  greater  depths  in  the  crust. 
Figure  8  shows  three  models  which  explain  most  of  the  fundamental-mode 
attenuation  coefficient  data  in  Figure  9.  There  is  little  difference  among  the 
predicted  attenuation  coefficients  predicted  by  the  three  models;  they  all 
predict  values  which  are  slightly  high  at  periods  of  8  and  10  s  and  slightly 
low  at  6  s,  but  are  in  good  agreement  with  all  other  data  points.  The  model 
which  yields  the  best  value  of  Qq  was  obtained  for  the  case  in  which  ^  =  0.5 
between  0.0  and  6.0  km.  That  model  also  predicts  a  value  of  y]  which  is 
within  one  standard  deviation  of  the  mean  observation.  When  $  =  1.0  at 
depths  between  0.0  and  4.0  km,  the  resulting  model  also  predicts  Qq  and 
T|  values  which  are  within  observed  standard  deviations.  Both  of  these 
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models  are  characterized  by  a  high-Q^  layer  overlying  a  lower  layer  at 
greater  depths.  The  third  model,  in  which  ^  =  0.2  at  depths  between  0.0  and 
8.0  km  predicts  a  reasonable  Qq  value  but  the  iq  value  is  somewhat  high. 

These  results  indicate  that  models  with  a  surface  layer  of  moderately 
high  values  overlying  a  low-Q^  layer  do  nearly  as  well  as  models  without 
such  a  layer  in  explaining  both  fundamental-mode  Rayleigh  wave  attenuation 
coefficient  data  and  Qq  and  iq  values  for  Lg.  In  an  attempt  to  decide  which  of 
these  two  types  of  models  is  correct  we  computed  synthetic  seismograms  for 
several  models  and  compared  them  to  observed  seismograms  which  were 
obtained  for  paths  entirely  within  the  Basin  and  Range  province  (Figure  10). 
Those  events  for  which  seismograms  were  available  are  listed  in  Table  1  and 
their  locations  are  shown  in  Figure  1.  As  shown  in  Figure  11,  the  model  in 
which  Qjj,  ~  60  in  the  uppermost  crust  and  Q^,  ~  1000  at  mid-crustal  depths 
produces  Lg  synthetics  for  which  the  envelopes  are  in  good  agreement  with 
the  observed  wave  forms  in  Figure  10.  The  only  significant  difference  is  that 
the  codas  in  the  observed  wave  forms  continue  to  greater  times  than  in  the 
synthetics.  Models  in  which  ~  3.0  x  10“^  (Q^,  ~  300)  in  the  upper  4  km 
of  the  crust  (Figtire  8)  produce  synthetic  seismograms  with  a  large  phase 
which  arrives  about  40  s  after  the  onset  of  Lg  at  a  distance  of  400  km  and 
about  90  s  after  the  onset  of  Lg  at  750  km.  No  phase  with  such  large  ampli¬ 
tudes  is  present  on  any  of  the  observed  seismograms  in  Figure  10.  The  same 
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late  phase  is  also  present  on  the  synthetic  seismogram  computed  for  a  dis¬ 
tance  of  400  km  for  the  model  in  which  ~  7.5  x  10“®  (Qj^  ~  130),  but  has 
attenuated  sufficiently  so  as  to  be  nearly  absent  at  a  distance  of  750  km.  In 
the  presence  of  noise  or  Lg  coda  produced  by  scattering  it  is  possible  that 
such  a  phase  could  be  present  but  be  undectable.  The  computations  indicate 
that  in  the  upper  crust  of  the  Basin  and  Range  cannot  be  as  high  as  300, 
but  could  be  as  high  as  130  or  slightly  more. 

The  models  plotted  in  Figures  4,  6,  and  8  all  give  values  at  1  Hz. 
Since  our  inversion  process  yielded  values  of  the  frequency  dependence  of 
Q"^,  we  can  also  determine  models  for  other  frequencies  where  intrinsic  Q 
dominates  seismic  wave  attenuation.  Figure  12  shows  Q^.  models  of  the 
Basin  and  Range  for  periods  of  1  s,  30  s,  and  100  s  for  the  case  when  g  =  0.1 
at  all  depths.  Q^,  becomes  lower  as  period  increases,  being  especially  low  (~ 
40)  in  the  upper  crust  at  a  period  of  100  s.  At  mid-crustal  depths,  where  Q^, 
is  very  high,  there  is  little  difference  between  the  models  for  the  three 
periods.  At  greater  depths  the  difference  increases  but  is  much  smaller  than 
that  at  shallow  depths. 

For  a  model  in  which  i  =  0.0  between  0  and  15  km,  and  t,  =  0.5  at  greater 
depths  (Figure  13),  the  models  are,  of  course,  identical  to  a  depth  of  15  km, 
but  show  a  large  decrease  of  with  increasing  period  at  lower  crustal  and 
upper  mantle  depths.  This  firequency-dependent  model  exhibits  the  kind  of 
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variation  we  would  expect  if  viscous  flow  occurs  in  the  lower  crust  and  upper 
mantle;  Q^,  values  would  continue  to  decrease  with  increasing  period  and  the 
anelastidty  observed  at  seismic  frequencies  would  be  expected  to  merge  into 
viscous  flow  at  very  large  time  scales. 

7  DISCUSSION 
Q  Variation  with  Depth 

The  Qp,  models  in  Figures  4,  6,  and  8  were  obtained  for  a  wide  range  of  ^ 
distributions  with  depth.  Most  of  those  distributions  produced  models 
which  are  similar  in  their  major  features  —  low  Q^.  values  in  the  upper  crust, 
a  high-Q^,  zone  at  mid-crustal  depths,  and  lower  Q^,  values  at  greater  depths. 
The  only  models  which  depart  from  that  distribution  are  those  in  which  5  is 
high  in  the  uppermost  crust.  Such  distributions  produce  a  layer  of  high  Qjj, 
values  in  the  upper  few  km  of  the  crust.  At  depths  of  8  km  and  greater, 
however,  those  models  closely  resemble  models  obtained  from  inversions 
using  other  I  distributions. 

Models  with  near-surface  layers  in  which  values  at  1  Hz  are  about 
200  or  more  can  be  ruled  out  on  the  basis  of  our  comparisons  between  syn¬ 
thetic  and  observed  Lg  seismograms.  Models  with  more  modest  near-surface 
high-Q^,  layers,  cannot,  however,  be  ruled  out  on  the  basis  of  either  our 
inversion  results  or  seismogram  comparisons.  If  we  assume  that 


-19- 


distributions  are  similar  throughout  the  western  United  States,  however,  a 
study  by  Hough  and  Anderson  (1988)  near  Anza  in  southern  California, 
seems  to  mitigate  against  such  a  layer.  They  found  Q  values  in  the  upper  5 
km  of  the  crust  (50-100)  to  be  substantially  lower  than  those  between  5  and  15 
km. 

Low  values  in  the  upper  crust  have  been  attributed  to  the  presence 
of  fluids  in  a  network  of  cracks  (MitcheU,  1975,  1980)  and  is  consistent  with 
both  laboratory  results  (Gordon  and  Davis,  1968)  and  theoretical  predictions 
(O'Coimell  and  Budiansky,  1977).  Our  results  imply  that  fluids  are  present 
to  10  km  or  more  in  the  crust.  The  high  values  at  greater  depths  in  our 
models  imply  that  those  fluids  are  absent  at  those  depths.  The  prediction 
that  those  fluids  are  absent  is  consistent  with  the  prediction  of  Bailey  (1990) 
that  free  fluids  have  a  short  residence  time  in  the  ductile  portion  of  the  crust. 

Q  Variation  with  Frequency 

The  models  which  seem  to  best  explain  all  available  information  are 
those  in  which  values  are  low  in  the  upper  crust,  increase  rapidly  to  high 
values  at  mid-crustal  depths,  and  decrease  at  greater  depths.  Inversions 
using  a  small  uniform  value  of  ^  or  inversions  in  which  I  is  small  in  the 
upper  crust  and  increases  to  relatively  large  values  at  mid-crustal  depth  pro¬ 
duce  such  models  (Figures  4  and  6).  The  former  £  distribution  prociuces 
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models  in  which  varies  modestly  with  frequency  at  all  depths  (or  not  at 
aU  for  5  =  0.0),  while  the  latter  produces  models  which  are  independent  of 
frequency  in  the  upper  crust  and  which  exhibit  a  large  variation  of  Q^,  with 
frequency  at  greater  depths.  The  lack  of  a  frequency  dependence  for  at 
upper  crustal  depths  is  consistent  with  laboratory  results  on  samples  which 
contain  a  network  of  fluid-filled  cracks  (Gordon  and  Davis,  1968)  and  is 
predicted  theoretically  for  models  in  which  fluid  can  flow  in  cracked  solids 
(O'Connell  and  Budiansky,  1977). 

At  lower  crustal  depths  at  100  s  is  much  lower  than  at  1  s,  being 
less  by  an  order  of  magnitude  at  a  depth  of  about  35  km.  We  would  expect 
that  type  of  behavior  to  occur  if  rock  in  the  lower  crust  undergoes  a  transition 
from  being  solid,  but  anelastic,  at  seismic  frequencies  to  exhibiting  viscous 
flow  at  large  time  scales.  This  model  predicts  that  seismic  waves  at  1  Hz  will 
propagate  efficiently  through  the  lower  crust,  but  that  long-period  waves  will 
attenuate  more  rapidly. 

It  is  interesting  to  compare  this  model  with  that  of  Toksoz  et  al.  (1990) 
obtained  using  Rg  waves  recorded  in  the  northeastern  United  States  and  both 
Rg  and  Lg  waves  in  Scandinavia.  Data  in  the  two  studies  were,  respectively 
at  frequencies  >  0.5  Hz  and  >  1.0  Hz;  their  model  therefore  pertains  to 
higher  frequencies  than  those  of  the  present  study.  Their  best-fitting  models 
require  (1)  low  Q^,  values  (<  100),  which  increase  with  frequency  as  f®  ^  in 
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the  upper  crust.  (2)  moderate  values  (100-500),  which  increase  with  fre¬ 
quency  as  ®  at  depths  between  2  and  10  km,  and  (3)  high  Q^,  values 
500),  with  indeterminate  frequency  dependence,  in  the  lower  crust.  The  pat¬ 
tern  of  increasing  Q^,  values  with  depth  is  similar  to  the  findings  of  the 
present  study,  as  well  as  earlier  studies  (e.g.  Mitchell,  1975,  1980);  the  distri¬ 
bution  of  frequency  dependence  is,  however,  drastically  different  from  that  of 
the  preferred  model  in  the  present  study.  The  differences  are  undoubtedly 
due  to  the  different  frequency  ranges  of  the  data  used  in  the  two  studies. 
Toksoz  et  al.  (1990)  explain  the  frequency  dependence  in  the  uppermost  crust 
as  being  produced  by  fluid  flow,  but  that  dependence  only  occurs  at  frequen¬ 
cies  above  a  critical  frequency  which  may  be  as  low  as  1  Hz  (Toksoz,  et  al., 
1987).  Since  our  work  utilized  data  at  frequencies  of  1  Hz  and  less,  that  fre¬ 
quency  dependence  is  not  expected  to  prevail  in  our  study.  The  f^  °  fre¬ 
quency  dependence  of  Q^,  at  depths  between  2  and  10  km  in  the  study  of 
Toksoz  et  al.  (1990)  is  explained  as  being  due  to  scattering.  Again,  at  the 
lower  frequencies  of  the  present  study  scattering  effects  are  expected  to  be 
much  less  important.  We  find  the  greatest  frequency  dependence  (in  our 
preferred  model)  to  occur  in  the  lower  crust,  a  region  where  the  data  of 
Toksoz  et  al.  (1990)  was  unable  to  provide  information. 


Relation  of  Q^,  to  Shear  Strength  in  the  Crust. 
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Smith  and  Bruhn  (1983)  modeled  the  rheological  behavior  of  the  crust  in 
the  Basin  and  Range  using  shear  stresses  calculated  at  various  depths  for  rock 
types  expected  to  occur  there,  using  thermal  parameters  from  Lachenbruch 
and  Sass  (1977)  and  assuming  a  strain  rate  of  10“^^  s“^.  They  obtained  a 
model  for  the  eastern  Basin  and  Range  which  included  three  brittle  layers, 
the  first  at  depths  between  0  and  8  km,  the  second  between  15  and  17  km, 
and  the  third  between  25  and  27  km;  ductile  material  lay  beneath  all  of  those 
layers.  Their  model  for  the  central  portion  of  the  Basin  and  Range  lacked  the 
deeper  two  brittle  layers.  Since  our  model  includes  only  one  rapid 
increase  with  depth,  it  appears  to  correspond  with  the  model  of  Smith  and 
Bruhn  (1983)  for  the  central  part  of  the  Basin  and  Range. 

The  transition  from  brittle  to  ductile  rock  in  the  Basin  and  Range  model 
of  Smith  and  Bruhn  (1983)  is  at  a  depth  of  8  km.  They  indicate  that,  because 
of  the  uncertainty  in  the  parameters  used  to  derive  their  models,  their  models 
should  be  considered  to  be  only  rough  guides  to  the  rheological  behavior  of 
the  crust.  In  another  study,  Jackson  and  McKenzie  (1983)  suggest  that  the 
brittle-ductile  transition  in  extensional  basins  occurs  somewhere  in  the  depth 
range  10-15  km  where  temperatures  are  about  350°C.  The  difference  between 
the  Smith  and  Bruhn  and  Jackson  and  McKenzie  models  can  be  explained,  at 
least  partly,  because  the  two  studies  make  different  assumptions  as  to  the 
temperature  at  which  the  brittle-ductile  transition  occurs.  Smith  and  Bruhn 
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take  that  temperature  to  be  less  than  300°C,  or  more  than  50°  cooler  than  that 
assumed  by  Jackson  and  McKenzie.  Sibson  (1984)  showed  that  several  fac¬ 
tors,  including  temperature,  austal  composition,  fluid  pressures  in  the  fric¬ 
tional  regime,  water  content  in  the  ductile  regime,  and  ductile  strain  rate, 
affect  the  depth  to  the  transition.  He  pointed  out,  for  instance,  that  the 
depth  to  the  brittle-ductile  transition  was  very  dependent  on  the  quartz  con¬ 
tent  of  rock  in  the  brittle  regime.  Smith  and  Bruhn  derived  their  model 
assuming  that  the  upper  crust  corresponds  to  Westerly  granite  which 
includes  about  30%  quartz  (Brace  et  al,  1965).  If  the  upper  crust  in  the  Basin 
and  Range  is  comprised  of  rock  with  a  lower  quartz  content,  the  brittle- 
ductile  transition  could  occur  at  depths  greater  than  8  km. 

Because  of  the  uncertainty  in  the  rheological  properties  in  the  Basin  and 
Range  we  have  compared  our  Q^,  with  two  models,  that  of  Smith  and  Bruhn 
(1983)  in  which  the  brittle-ductile  transition  occurs  at  a  depth  of  8  km  and 
that  of  Jackson  and  McKenzine  (1983)  where  it  occurs  in  the  depth  range  10- 
15  km. 

Figure  14  compares  the  model  of  our  study  with  the  model  of  shear 
strength,  (oj  -  03)72,  of  Smith  and  Bruhn  (1983)  and  with  a  model  for  an 
extensional  basin  in  which  the  brittle-ductile  transition  occurs  at  12.5  km,  the 
center  of  the  range  suggested  by  Jackson  and  McKenzie  (1983).  The  three 
models  are  aligned  with  a  histogram  of  earthquake  focal  depths  for  western 
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Nevada  as  compiled  by  Vetter  and  Ryall  (1983).  The  depths  are  taken  from 
Tables  2  and  3  of  their  paper. 

A  comparison  of  the  Q^,  model  with  the  shear  strength  of  Smith  and 
Bruhn  (1983)  would  suggest  that  only  the  range  of  decreasing  Q~^  values  (18 
to  13  X  10“^  between  about  3  km  and  8  km  depth  is  caused  by  either  crack 
closure  and  decrease  of  permeability  caused  by  increasing  pressure  (Brace, 
1980)  or  by  the  healing  of  cracks  at  elevated  temperature  (Smith  and  Evans, 
1984).  The  latter  process  may  occur  at  temperatures  as  low  as  200°C  in 
quartz.  At  depths  greater  than  8  km,  the  decrease  in  Q“^  (13  to  1  x  10“®) 
would  be  due  to  crack  closure  through  viscous  flow  as  well  as  to  decreased 
permeability  and  enhanced  healing  of  cracks  at  higher  pressures  and  tem¬ 
peratures.  If  the  brittle-ductile  transition  lies  at  greater  depths,  more  of  the 

decrease  is  caused  purely  by  increased  crack  closure  and  healing.  For 
the  example  in  Figure  14,  where  the  brittle-ductile  transition  is  at  a  depth  of 
12.5  km,  those  factors  cause  to  decrease  from  18  to  7  x  10“®  and  the  duc¬ 
tile  flow  contributes  to  the  Qj,,  decrease  from  about  7  x  10^  to  about  1  x  10“®. 

The  depth  distribution  of  earthquakes  in  the  western  Basin  and  Range 
(Figure  14)  appear  to  favor  a  rheological  model  in  which  the  brittle-ductile 
transition  is  in  the  range  10-15  km  as  suggested  by  Jackson  and  McKenzie 
(1983)  rather  than  at  shallower  depths.  Only  27%  of  the  earthquakes  occur  in 
the  upper  8  km  of  the  crust,  whereas  about  81%  occur  at  depths  of  14  km  and 
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less.  If  we  accept  that  model,  partial  crack  closure  by  pressure  and  annealing 
are  the  factors  which  produce  increases  in  Q^,  to  depths  of  10-15  km  and 
viscous  flow  is  an  additional  contributing  factor  at  greater  depths. 

At  depths  greater  than  about  25  km,  Q"'  values  increase  with  depth. 
The  effect  of  temperature  on  Q^,  at  those  depths  can  be  computed  using  a 
formulation  due  to  Anderson  and  Archambeau  (1964).  It  is  an  easy  matter  to 
describe  the  Q^,  decrease  using  that  formulation,  by  selecting  various  values 
for  activation  energies  and  activation  values  for  Q^,.  Correct  values  for  those 
parameters  are  not  known,  however,  so  we  will  not  present  the  results  of 
those  calculations,  but  simply  indicate  that  various  combinations  of  those 
parameters  can  achieve  the  desired  result.  Other  factors  which  may  contri¬ 
bute  to  the  increase  in  Q~^  at  these  depths  are  partial  melt  (Shankland  et  al, 
1981)  and  enhanced  dislocation  motion  (Minster  and  Anderson,  1980). 

8  CONCLUSIONS 

Various  assumed  distributions  of  frequency  dependence  produce  some¬ 
what  different  models  for  the  oust  of  the  Basin  and  Range  province. 
Many  models  explain  the  observed  attenuation  of  fundamentai-mode  Ray¬ 
leigh  waves  as  well  as  Q  and  its  frequency  dependence  for  1-Hz  Lg  waves. 
Most  models,  however,  have  the  common  features  of  low  Q^,  values  (50-80)  at 
depths  between  0  and  8  km,  rapidly  increasing  values  at  mid-crustal  depths 
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reaching  a  maximum  at  a  depth  of  about  20  km,  and  decreasing  values  at 
greater  depths.  Models  in  which  the  upper  crust  has  a  layer  of  higher 
values  (80-150)  overling  a  low-Q^.  region  can  also  explain  available  data,  but 
are  inconsistent  with  results  of  Q  studies  in  nearby  regions.  The  frequency 
dependence  of  Q^,  in  the  upper  crust  must  be  low  (0.0  -  0.1),  but  it  may  be 
substantially  higher  (as  much  as  1.0)  in  the  lower  crust.  A  model  in  which 
that  frequency  dependence  is  0.5  leads  to  a  model  in  which  can  be  an 
order  of  magnitude  lower  (~  100)  at  100  s  than  it  is  at  1  s  (~  1000). 

The  low  Q^,  values  in  the  uppermost  crust  are  most  easily  explained  as 
being  due  to  a  network  of  fluid-filled  cracks  which  permeates  brittle  rock  in 
that  depth  range.  Increasing  closure  of  those  cracks  by  pressure  and  anneal¬ 
ing  effects  causes  to  increase  with  increasing  depth.  They  are  the  dom¬ 
inant  effects  until  the  brittle-ductile  transition  is  reached  where  crack  closure 
due  to  viscous  flow  becomes  important.  The  depth  of  that  transition  is  not 
well  determined,  but  the  distribution  of  focal  depths  in  the  western  Basin  and 
Range  favors  its  occurrence  in  the  range  10-15  km.  Decreasing  Q^,  at  depths 
greater  than  25  km  may  be  caused  by  increasing  temperature,  increasing  pres¬ 
ence  of  partial  melt,  enhanced  dislocation  motion  or  some  combination  of 
these  effects. 
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FIGURE  CAPTIONS 

Figure  1.  Map  of  a  portion  of  the  Basin  and  Range  province  (denoted  by  the 
heavy  line)  showing  the  locations  of  stations  MNV  and  ELK  used  to 
obtain  fundamental-mode  Rayleigh-wave  attenuation  coefficients  used  to 
invert  for  internal  friction  models.  The  X  symbols  indicate  epicenters  of 
the  earthquakes  listed  in  Table  1  and  for  which  seismograms  are  plotted  in 
Figure  10. 

Figure  2.  Synthetic  seismograms  computed  by  modal  summation  (wang, 
1981)  for  distances  between  400  and  1200  km.  The  seismograms  were 
computed  along  an  azimuth  of  45°  from  a  vertical  strike-slip  fault  at  a 
depth  of  5  km  and  oriented  in  a  north-south  direction.  The  time  scale  is 
reduced  at  each  distance  R  using  a  reduction  velocity  of  5.0  km/s. 

Figure  3.  (Left)  model  at  1  Hz  obtained  from  the  inversion  of 

fundamental-mode  attenuation  coefficients  for  the  path  between  stations 
ELK  and  MNV  assuming  that  the  frequency-dependence  exponent  (Q 
is  0.0  at  all  depths.  Horizontal  lines  denote  one  standard  deviation. 
(Right)  Normalized  resolving  kernels  determined  for  depths  of  3.0,  9.0, 
19.5,  and  28.5  km. 

Figure  3.  S)nithetic  seismograms  computed  by  modal  summation  (Wang, 
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1981)  for  distances  between  400  and  1200  km.  The  seismograms  were 
computed  along  an  azimuth  of  45°  from  a  vertical  strike-slip  fault  at  a 
depth  of  5  km  and  oriented  in  a  north-south  direction.  The  time  scale  is 
reduced  at  each  distance  R  using  a  reduction  velocity  of  5.0  km/s. 

Figure  4.  Models  of  shear-wave  internal  friction  at  1  Hz  (left)  resulting  from 
the  inversion  of  fundamental-mode  Rayleigh-wave  attenuation  coefficients 
for  three  models  for  which  the  frequency  dependence  of  internal  friction 
(right)  is  uniform  with  depth  and  takes  on  values  of  0.0,  0.1,  or  0.2. 
Values  of  QLg  and  its  frequency  dependence  at  1  Hz  are  listed  for  each 
model. 

Figure  5.  Comparison  of  Rayleigh-wave  attenuation  coefficient  values, 
predicted  by  models  in  which  ^  is  uniform  with  depth  and  takes  on  values 
of  0.0,  0.5,  and  1.0,  with  observed  values. 

Figure  6.  Models  of  shear-wave  internal  friction  at  1  Hz  (left)  resulting  from 
the  inversion  of  fundamental-mode  Rayleigh-wave  attenuation  coefficients 
for  three  models  for  which  the  frequency  dependence  of  internal  friction 
(right)  increases  firom  0.0  to  a  positive  value  (0.2,  0.5,  or  1.0)  at  mid-crustal 
depths.  Values  of  QLg  and  its  frequency  dependence  at  1  Hz  are  listed  for 


each  model. 
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Figure  7.  Comparison  of  Rayleigh-wave  attenuation  coefficient  values 
predicted  by  the  model  in  Figure  6  with  observed  values.  The  solid  and 
dashed  lines  correspond  to  similar  lines  representing  the  models  in  Figure 
6. 

Figure  8.  Models  of  shear- wave  internal  friction  at  1  Hz  (left)  resulting  from 
the  inversion  of  fundamental-mode  Rayleigh-wave  attenuation  coefficients 
for  three  models  for  which  the  frequency  dependence  of  internal  friction 
(right)  decreases  from  a  positive  value  (0.2,  0.5,  or  1.0)  to  0.0.  Values  of 
QLg  and  its  frequency  dependence  at  1  Hz  are  listed  for  each  model. 

Figure  9.  Comparison  of  Rayleigh-wave  attenuation  coefficient  values 
predicted  by  the  model  in  Figure  8  with  observed  values.  Predicted 
values  for  the  three  models  are  nearly  indistinguishable. 

Figure  10.  Seismograms  from  the  broad-band  stations  ELK  and  MNV  for  the 
events  of  Table  1  and  Figure  1  after  being  high-pass  filtered  above  0.3  Hz. 

Figure  11.  Synthetic  Lg  seismograms  corresponding  to  paths  between  the 
events  shown  in  Figure  1  and  the  stations  ELK  and  MNV.  (Top)  Seismo¬ 
grams  predicted  by  the  model  in  Figure  6  in  which  Q~^  in  the  uppermost 
crust  is  about  17.5  x  10"^  (Q^^  ~  60)  with  higher  values  at  greater  depths. 
(Middle)  Seismograms  predicted  by  the  model  in  Figure  8  in  which  Q”^  in 
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the  uppermost  crust  is  about  3  x  10“®  (Q^,  ~  300)  with  lower  values  of  Q^,, 
at  greater  depths.  (Bottom)  Seismograms  predicted  by  the  model  in  Fig¬ 
ure  9  in  which  Q“^  in  the  uppermost  crust  is  about  7.5  x  10"^  (Qj,,  ~  130) 
with  lower  values  of  Q^,  at  greater  depths. 

Figure  12.  Frequency-dependent  model  of  the  crust  in  the  Basin  and 
Range  when  the  frequency-dependence  exponent  (Q  is  constant  with 
depth  and  has  a  value  of  0.1. 

Figure  13.  Frequency-dependent  model  of  the  crust  in  the  Basin  and 
Range  when  the  frequency-dependence  exponent  (£)  has  the  value  of  0.0 
at  depths  between  0  and  15  km  and  has  the  value  of  0.5  at  greater  depths. 

Figure  14.  a.  Histogram  of  earthquake  activity  at  2  km  intervals  in  the 
western  Basin  and  Range  (from  Vetter  and  Ryall,  1983).  b.  model  of 
the  Basin  and  Range,  c.  Model  of  shear  strength  obtained  by  Smith  and 
Bruhn  (1983)  for  the  central  part  of  the  Basin  and  Range,  d.  Model  of 
shear  strength  for  extensional  basins  by  Jackson  and  McKenzie  (1983)/  tak¬ 
ing  the  brittle-ductile  transition  to  be  at  12.5  km.  Stippling  represents  the 
brittle  portion  of  the  crust. 
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Lg  ATTENDATION  ACROSS  NORTHERN  EDRASIA 

We  have  mapped  Lg  coda  Q  and  its  frequency  dependence  throughout 
northern  Eurasia,  including  the  Barents  shelf.  Those  results,  when 
combined  with  a  similar  map  for  southern  Eurasia  (Xie  and  Mitchell,  1991) 
provide  Lg  coda  Q  and  frequency  dependence  values  which  cover  most  of 
the  continent,  from  the  Arctic  regions  of  Europe  to  southern  Asia  as  far 
south  as  20°  north  latitude.  Our  Q  map  for  Eurasia,  derived  from  data 
processed  up  to  now,  appears  in  Figure  1. 

Our  map  of  Eurasia  shows  Q  values  generally  decrease  from  north  to 
south  from  a  high  of  about  1000  to  about  200.  Low  values  (200-300)  occur 
in  a  broad  zone  which  stretches  across  western  China,  northern  India,  and 
adjacent  portions  of  central  Asia.  Localized  regions  of  low  Q  lie  within, 
and  near,  the  Caspian  Sea  and  western  Turkey.  The  highest  values  occur 
in  stable  regions  and  the  lowest  values  occur  in  regions  of  recent 
tectonic  activity,  results  which  are  consistent  with  our  earlier  studies 
in  other  regions. 

It  is  interesting  that  Q  values  for  the  region  of  the  Siberian 
platform  (higher  than  average  values  in  the  northeastern  portion  of  the 
map)  are  lower  than  those  for  the  East  European  platform  (high  values  in 
the  northwestern  portion  of  the  map) .  These  low  values  may  reflect 
intraplate  tectonism  in  that  platform  which  occurred  at  a  time  which 
postdates  the  accretion  of  the  subplates  which  formed  Eurasia. 

Work  is  continuing  on  this  map.  Data  from  a  recent  moderately  large 
earthquake  in  the  Netherlands  recorded  at  several  European  stations,  and 
data  from  Novaya  Zemlya  explosions  recorded  at  Arcess  and  Noress  are 
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currently  being  processed.  When  that  is  completed,  estimates  of  resolution 
and  error  in  the  measurements  will  be  obtained. 


